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The phase diagram of charge and color neutral two-flavor color superconducting quark matter is 
studied including the homogeneous two-flavor superconductor (2SC) and the inhomogeneous Fulde- 
Ferrell (FF) phases within the Nambu-Jona-Lasinio model. The low-temperature domain T < 
5 MeV of the phase diagram contains the FF phase, which borders at high temperatures to the 
2SC phase. The critical temperature of phase transition from the 2SC to the unpaired state is in 
the range 20-30 MeV. We derive the equation of state of matter and its composition and show that 
matter in mature compact stars should be in the inhomogeneous FF-like superconducting state. We 
briefly discuss the astrophysical implications of such a phase in compact stars. 
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dl' I. INTRODUCTION 
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I The quark substructure of baryons, predicted by quantum chromodynamics (QCD), suggests that nucleonic matter 
. may deconfine at sufficiently high densities attainable in compact stars. Compact stars featuring quark cores, known 
' as hybrid stars, may be distinguishable from their purely nucleonic/hyperonic counterparts by their astrophysical 
, manifestations, which may include magnetic, thermal, and rotational evolutions on secular timescales, as well as 
• shorter transients. 

'"^ , In this paper we derive the equation of state and composition of two-flavor dense quark matter at nonzero temper- 
^ ' ature under the conditions of charge neutrality and /3 equilibrium. At moderate densities quark matter is composed 
(— I of up and down quark flavors only, which in newborn protoneutron stars occur at near equal amount (the isospin 
'— chemical is vanishingly small). However, rapid compression of matter leads to its neutronization via the inverse beta 
' decay and the isospin chemical potential grows until it reaches the asymptotic value determined by the /3-equilibrium 
^ , condition Hd — fJ-u — Me among d and u quarks and electrons (here fii with i = d, u, e are their chemical potentials, 
J — ' respectively). Therefore, under stellar conditions the pairing between the two light flavors of quarks occurs at finite 
isospin chemical potential, i.e., when the Fermi surfaces of up and down quarks are shifted apart by an amount which 
T-H could be of the same order of magnitude as the gap in the quasiparticle spectrum. Under these conditions the actual 
pairing pattern may be significantly different from that of the BCS and is likely to involve breaking of the spatial 
\^ symmetries by the condensate order parameter, in analogy to less exotic low-temperature systems found in condensed 
matter (e.g. , electron gas or ultracold atomic vapor) 1-3]. 

The phase diagram of charge and color neutral quark matter, which features homogeneous condensates at non- 
zero temperature and moderate densities, is well established at the mean-field level [J-Q. It contains a sequence of 
^ , phases with the two-flavor superconducting (2SC) phase at low densities and the color-flavor-locked (CFL) phase at 
r asymptotically high densities (the color superconductingphases are reviewed in Ref . 0] ; for a review of astrophysical 
rS \ implications of color superconductivity see, e.g. , Ref. [3, Q). The domain occupied by the inhomogeneous phases 
. is difficult to access in general, since most of the studies are carried out in the regime where the Ginzburg-Landau 
■ " " ' expansion of the mean-field thermodynamic potential is valid [Tol - [T^ . An exception is the so-called Fulde-Ferrell 
(or plane- wave crystalline color superconducting) phase [l5l - [2l| . which admits straightforward treatment in the full 
domain of interest to compact stars. Low-dimensional treatments of inhomogeneous phases admit some simplifications 
as well ^25]. In effective four-fermion coupling models, like the Nambu-Jona-Lasinio (NJL) model, the interactions can 
be easily varied in a broad range of (effective) couplings that may account for medium modifications of the interactions 
in dense quark matter. 

In this paper we examine the phase diagram of charge and color neutral quark matter at non-zero temperatures and 
moderate densities and show that the simplest realization of the phases with broken spatial symmetries - the Fulde- 
Ferrell phase - arises in the low-temperature domain of the phase diagram. We do so, by extending the discussion 
in Ref. [2^ to include the conditions of charge and color neutrality, which impose additional constraints on the set 
of equations that are solved. As indicated above, we have chosen the spatial dependence of the order parameter 
that corresponds to a single plane wave. Therefore, our calculations offer only an upper bound on the ground state 
energy, since more complicated order parameters have been proven to possess lower energy in the domain where the 
Ginzburg-Landau expansion is valid. However, qualitatively the phases with broken spatial symmetries share many 
common features and the Fulde-Ferrell (FF) phase is a convenient starting point of exploration. 

Our discussion is restricted to the two-flavor case appropriate at moderate densities, i.e., its validity domain is 
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constrained by the temperature and density where strange quarks appear in substantial amounts. The role of strange 
quarks is twofold: first, because of their negative charge it becomes energetically favorable to replace high momentum 
electrons by low momentum strange quarks; second, strange quarks disfavor the BCS type cross-flavor pairing among 
the light quarks and favors instead the FF or other crystalline-type pairing. If the Fermi energy of strange quarks 
is large enough they may choose to pair in a color-flavor-locked pattern which, except for the asymptotically high 
densities, will involve some mismatch between the Fermi surfaces with possibly multiple FF phases involving quarks 
of all three flavors. 

Finally, a number of non-crystalline phases have been suggested in the literature, which address the stress caused 
by the mismatch of Fermi surfaces differently. These alternative phases invoke separation into normal and supercon- 
ducting domains jHI or pairing induced changes in the shapes of the Fermi surfaces [2^. 

This paper is organized as follows. Section HIl presents the key equations describing the FF phase (see also Ref. 20\). 
In Sec. mil we discuss the conditions of charge and color neutrality. Our results are presented in Sec. IIVI and our 
conclusions are collected in Sec. |Vl 



II. GAP EQUATIONS AND THERMODYNAMICS OF FF PHASE 



The mean field (MF) thermodynamic potential of the FF phase is given by 

JImF = i^MF + ^MF) (1) 

where the first term includes the contribution from the condensate of red-green quarks and the second term is the 
contribution from the unpaired blue quarks. The first term is given by 

a2a2 f H^k f -1 

^MF = 3^-2^ J j^{Et{A)+E-{A)-2T\nf[~Et{A)]-2T\nf[~E-{A)]}, (2) 

where A is the cutoff of the NJL model, g is the strong coupling constant, /(w) = [1 +exp(a;/T)] ^ is the Fermi 
distribution, and T is the temperature. The contribution of the blue quarks is 




where / summation is over flavors and /i^ is the chemical potential of blue quarks of flavor /. The thermodynamic 
potential depends on the quasiparticle energies, which are given by 

EtiA) « Ea,,±^EI^ + \A\^, (4) 

where Es.e and EA,e are the parts of the spectrum which are even (symmetric) and odd (asymmetric) under exchange 
of quark flavors in a Cooper pair. These are given by^ 



EsA\k\AQ\,0,p) = E+-efl, (5) 
EA,emAQ\,d,Sti) = S^i + eE-, (6) 

with E^ — {\k + Q/2\ ± |A; — Q/2|)/2, where 6 is the angle formed by vectors k and Q, Sfi = (/i„ — fid)/2 and 
A = (Mm + Md)/2. The spectrum of unpaired blue quarks is given by ^^(fe, fJ,b) = Es^k, 0, 0, jib) ~ \k\ — enb, where ^b 
is the chemical potential of blue quarks. 

The form of quasiparticle spectrum requires a subtraction from the thermodynamic potential of the unphysical 
contribution from currents that exist in the case where pairing gap vanishes. Therefore, we need to evaluate 

SImf — f^MF — f^MF° + ^mf"'*^~"- C^) 



Because \Q\ Sfi is small compared with the typical scale fi for quark momenta, we can keep the leading order in \Q\/fJ- terms, in which 
case the spectrum takes the form given in Ref. j20j. 
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The stationary point(s) of the thermodynamic potential determine the equihbrium values of the order parameters: 



MF 



OA 



-0, 



dn 



MF 



d\Q\ 



0. 



(8) 



The direction of the vector Q, i.e., the axis of symmetry breaking is chosen by the system spontaneously; the 
appearance of a preferred direction implies that 0(3) symmetry of the BCS theory is broken down to 0{2) in the FF 
phases; phases with more complicated lattice-type order parameters break the translational symmetry of the problem. 
The relative concentrations of the u and d quarks are determined from the derivative of the thermodynamic potential 
with respect to the corresponding chemical potential as 



dfl 



MF 



MF 



OfJ-u OfM 

The explicit form of the gap equation, which follows from Eqs. ([2]) and 



IS 



tanh 



tanh 



(9) 



(10) 



Note that in the isospin-symmetric limit the quasiparticle spectra are degenerate ii'g'"(A) — -~E^{A) and (with 
IQI = 0) the gap equation reduces to the ordinary BCS gap equation for ultrarelativistic fcrmions. Finally, the 
equation for the total momentum is given by 



d^k (dEti^ 
-, — 7^ \ — ^, „, — tann 
(27r)3 \ d\Q\ 



2T 



dE-{A) 
d\Q\ 



tanh 



E7{A) 



2T 



(11) 



with 



d\Q\ 



e ^ Es,em,\Q\,0,f^) 
2 E+{A)~Ee{A) 



Q + 2k cos 9 
2\2k + Q\ 



e EsM\kUQ\Afi)_ 
2^ E+{A)~Ee{A) 



Q - 2k cose 
2\2k~Q\ 



(12) 



The zero temperature counterparts of Eqs. ([TU|) -([TT |) are straightforward to obtain with the help of the identity 
tanh(x/2) = 1 — 2f(x) and the limiting form of the Fermi distribution function /(x) — 6{~x) for T — > 0, where 9{x) 
is the Heaviside step function. 

As long as strange quarks are too heavy to appear in matter, the net positive charge of quarks is neutralized 
by electrons. At relevant densities and temperatures the electrons can be treated to a good approximation as a 
noninteracting gas. Their thermodynamic potential is then given by 



^cicc — 2T 



[In/ (-Gec(fc.Me)) +ln/(-C,,(fc,A*e))] 



(13) 



with £,^if,^{k, fie) = y/k'^ + ml — /Xe, where me is the electron mass. The second term in Eq. represents the 
contribution of anti-particles (positrons) to the thermodynamical potential with C(Jcc(^' Me) = \/k'^ + m^ + fie- 



III. IMPOSING CHARGE AND COLOR NEUTRALITY 



The matrix representing the chemical potentials of quarks in the color-flavor space can be written in terms of the 
color group generators and the matrix of electrical charge operators as 



where /is is the baryon chemical potential, the electrical charge operator in flavor space is 

'2 1 1^ 



Qe = diag 



3' 3' 3/' 



while the color operators are given by 



= idiag(l,-l,0), = ^diag(l,l,-2). 



(14) 



(15) 



(16) 



4 



We choose the direction of the color symmetry breaking the blue direction. The densities of red and green quarks are 
then automatically equal, therefore = and there is no need for the conjugate to chemical potential ^3. From 
Eqs. (Il4])-(I16|), 



Ms _ Me 
3 6 

Me 



Ms 

2V3' 



(17) 
(18) 



The full thermodynamic potential of our system is now given hy fl — fluF + ^^cicc- By choosing instead of chemical 
potentials of quarks the chemical potentials /xg s-^d Me as independent thermodynamic variables we obtain a set of 
equations that should be solved simultaneously: 



on 



MF 



dA 



= 0, 



on 



MF 



d\Q\ 



= 0, 



(19) 



d^JLe 



0, 



0. 



(20) 



The densities of quarks and leptons can be derived through derivatives of the thermodynamic potential with respect 
to the conjugate chemical potentials; for electron density we have 



dVL. 



nc\c 



elec 



= 2 



dUe ~J (2^)3 

The densities of paired up and down quarks are given by 

d^k eEs.e 



[/(Cec)-/(Gcc)]- 



rq I ro 



-E 

e 

E 



tanh 



2T 



tanh 



2T 



d^k 
(2^ 



tanh 



^e+(A) 

2T 



tanh 



E7{A) 
2T 



The densities of unpaired up and down quarks are given by 



2E 



e/(6(fe,M^)), 



(27r)3 



The total densities of up and down quarks are 



The baryon density is 



nd 



rq , b 



n,, = ri„ 



riB 



riu + Ud 



(21) 

(22) 
(23) 

(24) 
(25) 



(26) 
(27) 



(28) 



The electrical and charge neutrality require, respectively. 



2 1 

- - "dec = 0, 



0. 



n7 + rC-2n^-2r4 



(29) 
(30) 



The thermodynamic properties of matter are completely determined through four gap equations (|19p -(|20 p for a given 
temperature and baryon density (or baryon chemical potential). 
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TABLE I: The column entries are as follows: baryonic density in units of fm~ , total pressure and electron pressure in units of 
MeV fm"'', the baryonic chemical potential, the chemical potentials of u, d quarks and electrons in units of MeV, the densities 
of d quarks, u quarks, and electrons in units of fm~^, the gap in units of MeV, and total momentum of the condensate in units 
of fm-i. 
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0.29 


0.40 


90.2 


37.6 


1422.0 


539.4 


512.9 26.5 
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r 


= 5 MeV 












0.32 


66.5 


29.0 


1322.5 


501.5 


475.8 25.7 


0.44 


0.55 


0.14 


37.8 


0.00 


0.36 


77.1 


32.8 


1369.4 


519.2 


493.3 25.9 


0.49 


0.61 


0.16 


45.7 


0.00 


0.40 


90.7 


37.6 


1420.1 


538.6 


512.0 26.6 


0.55 


0.67 


0.18 


53.7 


0.00 


0.44 


105.6 


42.7 


1468.4 


557.2 


529.7 27.5 


0.61 


0.74 


0.20 


61.8 


0.18 


0.48 


121.7 


48.0 


1514.7 


575.1 


546.8 28.3 


0.67 


0.80 


0.21 


70.2 


0.29 



T = 10 MeV 



0.32 62.7 26.7 1309.5 495.6 474.6 21.0 0.43 0.53 0.14 35.9 0.00 

0.36 86.6 39.3 1370.8 523.1 488.0 35.2 0.46 0.62 0.16 43.6 0.00 

0.40 100.4 44.1 1422.8 542.5 507.7 34.8 0.53 0.69 0.18 52.7 0.00 

0.44 111.3 46.3 1468.9 558.8 527.6 31.1 0.60 0.74 0.19 61.4 0.00 

0.48 126.1 50.5 1513.9 575.6 545.2 30.4 0.66 0.80 0.21 70.0 0.00 



IV. RESULTS 

Eqs. ([T9|) and ([20| were solved numerically by finding a solution that satisfies simultaneously the four integral 
equations. The ultraviolet divergence of the integrals was regularized by a three-dimensional cut-off in momentum 
space \p\ < A. The phenomenological value of the coupling constant G = jVlts? in the (gg) Cooper channel is 
related to the coupling constant Gd in the (gg) di-quark channel by the relation G — Nc/{2Nc — '2)Gd, where Nc = 3 
is the number of quark colors. The coupling constant Gd and the cut-off A are fixed by adjusting the model to the 
vacuum mass and decay constant of the pion. As in Ref. [1^ we employ the parameter values Gd = 3.11 GeV~^ and 
GdA2 = 1.31. 

The phase diagram of charge and color neutral matter is shown in Fig. [T] The baryonic density range extends up to 
3poi where po = 0.16 fm~"^ is the nuclear saturation density. The model is limited to the above density range, because 
at somewhat larger densities the Fermi momentum becomes comparable to the cut-off of the model. The covered 
temperature range extends up to the critical temperature for the onset of color superconducting phase transition. 
It is seen that the low-temperature domain is occupied by the FF phase; in the intermediate temperature regime 
5 < T < 20 MeV, the FF phase is replaced by a homogeneous 2SC phase, the phase transition being first order. At 
the critical temperature, which depending on the density lies in the range Tc ~ 20 — 30 MeV, there is a second order 
phase transition from the 2SC to the unpaired phase. 

The equation of state of matter and related parameters are tabulated in Table HI The large pairing gaps at large 
densities are due to the density scaling of the density of states at the Fermi surface; note that the magnitude of 
the gap depends on the chosen value of the cutoff and is, therefore, specific to the present model. Furthermore, the 
pressure requires some normalization procedure, which is chosen to be the standard NJL prescription, which requires 
zero pressure (and energy density) in vacuum (i.e., at T = /j, = 0). In practice, the low density limit of the model 
should correspond to quarks confined in baryons, so that the normalization should be controlled by the confinement 
mechanism in non-zero temperature/density QCD. This freedom allows one to add a bag constant to the pressure 
of the NJL model to match it to the equation of state of hadronic matter at a desired density and temperature 
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FIG. 1: (color online). The phase diagram of neutral two-flavor quark matter. The high-temperature low-density region 
corresponds to unpaired quark matter, the intermediate temperature region corresponds to the homogeneous 2SC phase, and 
the low-temperature domain corresponds to the inhomogeneous FF phase of quark matter. 



of deconfinement. The point of onset of deconfinement is, thus, a free parameter of the model. Consequently, the 
pressure of the model may be shifted by a constant value of the bag. 

As seen from Tabled the equation of state is only weakly temperature dependent. Table |T] shows also the composition 
of matter and the chemical potentials of species, which are related by the /^-equilibrium condition /i^ = /ie + Mm- 

This 

condition will be violated in the high temperature domain T > 10 MeV of our phase diagram, where neutrinos are 
thermally populated as a consequence of their short mean-free path at such temperatures. The non-zero neutrino 
chemical potential then must be added to the left-hand side of the /3-equilibrium condition (|18p . Since there is no or 
little overlap between the domain occupied by the FF phase and the domain where neutrinos are trapped, we shall 
not repeat our analysis at non-zero neutrino chemical potential. Let us note that the phase transition from the 2SC 
to the FF phase is first order and is associated with a latent heat; such a heat release in compact stars during the 
collapse stage may serve as an engine that powers high-energy burst phenomena from the collapsing star. 



V. SUMMARY 



In this paper we have investigated the phase diagram of charge and color neutral quark matter in two-flavor NJL 
model. Our results strongly suggest that the 2SC phase does not exist in the interiors of mature compact stars. 
Instead, the stellar matter at intermediate densities is in the FF, or related crystalline, phase. A possible alternative, 
which is not excluded by our calculations, is that nuclear matter directly transforms to three-flavor CFL matter. 
The FF phase provides an upper bound on the energy of the color superconducting phase, which can presumably be 
lowered by choosing more complicated forms of the order parameter. The model parameters in our study were chosen 
to correspond to the common values adopted in the NJL model; in the unlikely case where the in-medium couplings 
are substantially larger, the 2SC phase may have lower energy than the FF phase (see Ref. [20j). 

Provided a nuclear equation of state and a matching procedure {e.g. , Maxwell construction) the present equation 
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of state can be used to construct the coniplete equation of state of matter in compact stars and their models (such 
constructions are reviewed in, e.g. , Refs. [H, H^])- The interesting possibihty that three-flavor crystalline phase may 
appear in stable massive compact stars has been pointed out recently [1^. The presence of a crystallinelike phase in 
the interiors of hybrid compact stars may have a number of interesting astrophysical implications; e.g. , the finite shear 
modulus of matter implies deformations of the quark core that can generate gravitational waves at the level of the 
sensitivity of current gravitational wave observatories [l!] • The properties of the FF phase may manifest themselves in 
the cooling of compact stars. The emission rate of the dominant Urea process has been computed in phases similar to 
the FF phase in Refs. [29, 3C|. Finally, the first order BCS-FF phase transition and the latent heat release associated 
with it may power high-energy burst phenomena from collapsing stars. 
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